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Abstract 

Analogy between an approximate version of Feigenbaum renormalization group analysis in 
complex domain and the phase transition theory of Yang-Lee (based on consideration of formally 
complexified thermodynamic values) is discussed. It is shown that the Julia sets of the renormal- 
ization transformation correspond to the approximation of Mandelbrot set of the original map. 
New aspects of analogy between the theory of dynamical systems and the phase transition theory 
are uncovered. 
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1 Introduction 

Transition between magnetic and non-magnetic 
phase of matter is intensively studied for a long 
time [U El El ll| • At low temperatures, elementary 
atomic magnets tend to be parallel providing co- 
herent (magnetic) phase of the system. With in- 
crease of the temperature, thermal fluctuations 
break this order drastically. If the temperature 
is higher than the phase transition critical point 
(Curie point), then chaotic orientation of elemen- 
tary magnets takes place (non-magnetic phase). 
However, near the critical point, the elementary 
magnets still keep order on certain range of dis- 
tances and over certain time intervals, which de- 
crease with growth of the temperature. Thus, at 
the critical point, areas of coherence and areas 
of thermal fluctuations coexist. Areas of fluctua- 
tions in the magnetic phase or areas of coherence 
in the non-magnetic phase may be distingueshed 
only within a certain range of scales. For large 
enough scales the system looks completely or- 
dered (as at zero temperature) or chaotic (as at 
infinite temperature). If one considers the same 
system at a definite temperature but within var- 
ious scales, it looks as the system at different 
temperatures. The transformation of scales cor- 
responds to the temperature renormalization. 

Let we have a cubic lattice of N atoms with 
interatomic distance a at temperature T. If one 
considers this system at rough enough scale, at 
which the elementary cell has the characteristic 
size na and contains n 3 atoms, the system looks 
as a lattice of N/n 3 atoms at the renormalized 
temperature T 1 = R(T) (R is a function of renor- 
malization). At the moment of the phase tran- 
sition, one can observe the thermal fluctuations 
within any scales. These fluctuations are self- 
similar. The idea of self-similarity has been of- 
fered by Kadanoff |5| and further advanced by 
Wilson to the renormalization method jj| (see 
also 0). 

Concepts of universality and scaling near the 
critical point and a renormalization group (RG) 
method have been transferred to nonlinear dy- 
namics by Feigenbaum IB]- The analogy be- 
tween transition to chaos and phase transitions 
consist in increase of temporal scales with ap- 
proach to a critical point, self-similarity, an op- 
portunity of application of RG analysis, and in 



existence of universal critical indexes and scaling 
factors, which are determined by the most general 
requirements to the type of system EH HU E2| • 

In 1952 Lee and Yang advanced an approach 
in the phase transition theory based on consider- 
ation of analytical properties of some thermody- 
namic values, such as partition function and free 
energy, depending on the temperature considered 
formally complex variable [T8| ITi] . 

The partition function is defined as follows: 

Z JV (T,F) = Eexp (-^) = 

(j \ (1) 
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where E is a total energy of a configuration of 
the system, H is an external magnetic field, T is 
a temperature, k is the Boltzmann constant, s« - 
is the spin variable, taking place at z-th cell of the 
lattice. Spins interact via magnetic field with the 
nearest neighbours (J > corresponds to ferro- 
magnetic interaction, J < - antiferromagnetic). 
Partition function can be transformed to 

N 

Z N = e- N / 2 J2Pn* n , (2) 

n=0 

by the variable change z = exp(—2H/kT). It 
is known, that at the phase transition the diver- 
gence of the thermodynamic potential (free en- 
ergy) 

kT 

f(H,T)= lim -^\nZ N (3) 

N^oo iv 

takes place. 

Obviously, zeros of the partition function (roots 
of the algebraic equation Zn = 0) are the candi- 
dates for the phase transition points. However, 
these roots should be positive and real, which is 
impossible as all factors p n > 0. In Yang-Lee the- 
ory, the temperature and the magnetic field are 
considered as complex variables [U 113 EE3| • Com- 
plex zeros of the partition function (called the 
Yang-Lee zeros), gather near the real axis. It be- 
comes possible only in the thermodynamic limit, 
the asymptotics of infinite number of atoms. 
When the number of atoms is finite, there is a fi- 
nite number of zeros in the complex plane. With 
increase of number of atoms, the set of zeros be- 
comes dense and nestles to the real axis more and 
more. In the thermodynamic limit the Yang-Lee 
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zeroes form a fractal set, that cross the real tem- 
perature axis at the point of the phase transi- 
tion. To find the Yang-Lee zeros it is appropri- 
ate to apply the Wilson renormalization method, 
which corresponds to a consecutive decrease of 
the number of degrees of freedom of the parti- 
tion function. It is necessary to find the trans- 
formation of the iV'-particle partition function to 
the iV-particle one (N' < N). This transforma- 
tion is not reversible. Thus, to get zeros of Zn 
it is necessary to have zeros of Zw, and then to 
construct their backward images using the RG 
transformation. Having repeated this procedure 
many times one comes finally to the trivial two- 
atom partition function. If zeros of Z2 belong to 
a basin of attraction of a fixed point of the RG 
transformation, then zeros of Z^ in the thermo- 
dynamic limit will coincide to the boundary of 
the basin of attraction in the complex plane. As 
shown by Derrida [T5] . this boundary, which is 
the Julia set of the RG transformation, is identi- 
cal to the set of Yang-Lee zeros. Unfortunately, 
the renormalization transformation may be per- 
formed analytically only for simple class of mod- 
els - the hierarchical lattices (for example, one 
and two-dimensional model with Ising spins 
hierarchical lattices with Potts spins [T7] I18j). 

Yang-Lee theory appears to be fruitful for un- 
derstanding the phase transitions. It seems that 
a similar approach to the analysis of transition 
to chaos in dynamical systems would be useful 
for deeper understanding of the analogy with the 
phase transitions, and for development of the new 
criteria describing complexity of behavior of non- 
linear systems. In the present paper we develop 
such approach basing on the approximate RG 
analysis of transition to chaos through the period 
multiplication bifurcations. 

2 Approximate RG analysis 

In 1978 Feigenbaum discovered universality of 
cascade of the period-doubling bifurcations and 
described it on a basis of the RG method. The 
simplest example representing the Feigenbaum 
universality class is quadratic map 

Xn+l = f\(Xn) = A - xl, (4) 

where a; is a real dynamical variable, and A is a 
real parameter. This map has fixed points, which 



can be found as roots of the equation = A — x 2 . 
If A > A x = 3/4, a cycle of period 1 (that is the 
fixed point) loses its stability. This Ai is the pa- 
rameter, at which the first period-doubling bifur- 
cation occurs (the multiplier of the fixed point is 
/1 = f\(x*) = —2x* = —1). Values of A for se- 
quent bifurcations can be found by means of ap- 
proximate RG method Let's apply the map 
(4) two times: 

x n+2 = A - A 2 + 2\xl - x 4 n (5) 

and neglect the last term, the fourth power of x n . 
Then, by the scale transformation 

Xn->Xn/a , a = -2A. (6) 

this map can be rewritten in the form x n+ 2 = Ai- 
re 2 , which differs from (4) only by renormalization 
of A 

A x = ip{\) = -2A(A-A 2 ). (7) 

Thus, the operator of evolution for the double 
interval of discrete time can be reduced to the 
original operator by the renormalization transfor- 
mation (7). Repeating this procedure with scale 
factors a% = — 2Ai, one can obtain a sequence 
of the same form 

%n+2 m — A m — x 2 , A m = </?(A m _i). (8) 

Fixed points of these maps correspond to the 
2 m -cycles of the original map (m = 1, 2, 3, ...). It 
is easy to see, that all these cycles, as well as the 
fixed point of the map (4), become unstable at 
A m = Ai = 3/4. Solving a chain of the equations 

Ai = ^(A 2 ), A 2 = <p(A 3 ), A m _! = ip(A m ) (9) 

we get the corresponding sequence of bifurcation 
values of parameter A (with A ~ A m the 2 m -cycle 
of (4) arises). From iteration diagram of Fig.l 
it is evident, that this sequence converges with 
m — > 00 to a definite limit Aoo, the fixed point of 
the RG transformation. It satisfies the equation 
Aoo = ^(Aoo), thus Aoo = (l + v^/2 « 1.37. The 
scaling factors also converge to the limit: a m — > 
a, where a = — IK^ « 2.74. The multipliers 
(Floquet eigenvalues of the 2 m -cycles) converge 
to the universal value \i m — > /i — \fl — 4Aoo ~ 
-1.54. 
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Table 1: Rigorous and approximate values of 
Feigenbaum critical indexes 





rugoious 


A nrvvAViTTi o i o 

r\ppi UXlIIlcLLt; 


A* 


1.401 


1.37 


a 


-2.802 


-2.74 


S 


4.669 


5.73 




-1.569 


-1.54 



From transformation (8) it is possible to ob- 
tain the law of convergence of the bifurcation se- 
quence: 



A m = v?(Aoo) + y?'(A.oo)(A m +i - Aoo) = 

= Aoo + 5(A m+1 - Aoo) 

(10) 

where 5 = (p'(Aoo) = 4 + v3 ~ 5.73 is a con- 
stant, characterizing the convergence to the crit- 
ical point. 

In table 1 we summarize the values of criti- 
cal indexes (critical point, scale factor, parame- 
ter scaling constant and universal multiplier) ob- 
tained by means of the exact and approximate 
RG analysis. The correspondence between them 
is well enough. 
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Figure 1: Iteration diagram of the RG transfor- 
mation (7). Dashed line designates the backward 
iterations starting at the first period doubling bi- 
furcation point (Ai = 3/4) and mapping to the 
further bifurcation points A m . 



3 Complex RG transforma- 
tion and its Julia set: 
The critical phase transi- 
tion line 

Let us regard the variable x and parameter A of 
map (4) as complex numbers. Then it becomes 
clear, that Feigenbaum' s universal scaling laws 
correspond to a special case of general scaling 
properties of Mandelbrot sets. The Mandelbrot 
set shown in Fig.2, can be determined as [I] 

M = {A G C : lim /™(0) + oo}. (11) 

n— »oo 

Thus, M is the set of values of complex parame- 
ter, for which the trajectories launched from the 
extremum of the map (4) — > A — > A — A 2 — , , 
remain in a finite domain. Mandelbrot set in- 
cludes parameter values corresponding to peri- 
odic trajectories of the map (the main cardioid 
and a number of the round "leaves", painted by 
gray color), and a set of parameter values cor- 
responding to bounded chaotic dynamics (black 
fractal pattern). Feigenbaum's cascade of period- 
doublings takes place along the real axis (see the 
picture of the bifurcation tree in the Fig.2). 

However, in the complex plane one can observe 
many other accumulation points of other bifur- 
cation cascades. One of such points, with scal- 
ing properies distinct from the Feigenbaum laws 
and intrinsic exclusively for complex analytical 
maps [2HII2I] is associated with the cascade of pe- 
riod tripling (221 • Also, there exist critical points 
connected with cascades of period quadrupling, 
period 5-tupling, etc. 

Thus, in the complex domain, a number of new 
scenarios of transition to chaos and other critical 
phenomena occurs. Let us consider the approx- 
imate RG transformation generalized on a com- 
plex plane 

A„+i = (f(X n ) = -2A n (A„ - \ 2 n ). (12) 

Starting from any point of a complex plane the 
trajectory of the map approaches the attractor at 
A = or escapes to the attractor at infinity. The 
borderline of these basins of attraction is the Ju- 
lia set J of the of the complex mapping (12) jl]. 
Let us call it the critical set by analogy with the 
theory of phase transitions, where it is interpreted 
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Figure 2: Correspondence between the Mandel- 
brot set (a) and the Feigenbaum bifurcation tree 
(b). Mandelbrot set is built up on a plane 
(ReA, ImA) for the complex quadratic map (4). 
Bifurcation tree is shown in the plane (A, x n ) for 
the map (4). Gray color in panel (a) indicates 
areas of periodic dynamics (the periods of cycles 
are marked); black color designates points cor- 
responding bounded in the phase space chaotic 
dynamics; white color means escape to infinity. 

as the critical phase border. Julia set J is deter- 
mined as the border of the set 

P = {A e C : lim <^ n (A) ^ oo}. (13) 

n^oo 

In Fig. 3 a diagram ot the complex plane of ini- 
tial values of A of RG transformation (12) is rep- 
resented. Black color indicates the basin of at- 
traction P, that is the area bounded by the Julia 
set J. Shades of gray color designate areas of dif- 
ferent time of escape (i.e. dynamic distance up to 
attractor at infinity), i.e. the sets of points with 
various number of iterations necessary to escape 
out of a large enough ball. 

It is evident that the Julia set in Fig. 3 is quite 
similar to the Mandelbrot set of the map (4) 
(Fig. 2). Similarity between them can be ex- 
plained as follows. Indeed, the escape of A to in- 
finity with iterations of RG transformation means 
the escape to infinity of iterations of the map 
(4). Distinctions can be explaned by approxi- 
mate character of RG transformation. Also, it is 
necessary to take into account that transforma- 
tion (12) describes properties only for the period- 



doubling cascade of period-doublings, thus pro- 
viding the similarity of sets M and J only near 
the Feigenbaum point. (The asymptotic simi- 
larity of the area of an analyticity of exact so- 
lution of the RG equations and Mandelbrot set 
near the Feigenbaum point was discussed earlier 
in works [2HJ I2U •) It is interesting to gener- 
alize these results for other sequences of period- 
multiplication bifurcation cascades. 




Figure 3: Complex plane of initial values of A 
for RG transformation (12) is represented. Black 
color designates the basin of attraction of origin 
bounded by the Julia set. Shades of gray color 
mark the area escape to infinity. Light colors cor- 
respond to the larger dynamical distance from the 
attractor at infinity. 

4 Complex cycles of RG 
transformation 

Let us describe procedure of visualisation of crit- 
ical set J with the help of a backward iterations. 
This methode is based on known fundamental 
properties of the Julia sets of analytical maps 

If there is a repeller point A* for the map A n+ i = 
<p(\ n ), then the set 

J' = {A: vAA) = A„n = l,2,...}, (14) 

is dense in J. The fixed point can be found di- 
rectly from an iterative polynomial, root of 
the equation (p(\*) = A*. If the fixed point satis- 
fies a condition |</?'(A*)| > 1, then it is repeller. 

Julia set is invariant both in respect to for- 
ward and backward iterations. Backward orbit 
of any point laying in some basin of attraction 
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approaches closer to a border of this basin, i.e. to 
the Julia set. 

Repeller points in our case are known to be 
A* = (1 ± y / 3)/2. The backward transformation 
also can be found analytically by solving equation 
(12) in respect to the variable A n (we easily do it 
with the help of software packet Mathematica) 



y = I + 22/3 | 

n 3 3(4+27A n+ i+3 v /3A n+ i(8+27 A n+ i)) 1 / 3 
(4+27A n+1 +3 v /3A n+1 (8+27A 

3x2 2 /3 ) 

^2,3 _ 1 l±V3i 

n 3 3x2 2 / 3 (4+27A n+ i+3 v /3A n+ i(8+27A n ~) 1 / 3 
(l=fcy / 3i)(4+27A n+ i+3 v /3A n+ i(8+27A n+ i)) 1 / 3 
6x2 2 / 3 

(15) 

Also the Julia set J can be regarded as a set 
of all unstable cycles of every possible periods 
n jU EB1 E2j. With increase of n this set be- 
comes more and more dense, and at n — > oo the 
distribution of the points on the complex plane 
yields the Julia set. Numerical calculation of 
an unstable cycle of period n is connected with 
constructing of every possible periodic sequences 
ei, E2, s n , E\, 62, s n , where £j takes values 
1, 2 or 3, which correspond to a choice of a root 
A 1 , A 2 or A 3 at the i-th backward iteration of (15). 
As a result of iteration procedure, the sequence of 
values of A converges to a cycle of period not ex- 
ceeding n. 

So, one of the unstable fixed points of RG trans- 
formation (12) is situated on the real axis and cor- 
responds to the Feigenbaum transition to chaos. 
Beside it, there are unstable cycles of any other 
period. Elements of a cycle of period n of RG 
transformation (12) correspond to the accumula- 
tion points of bifurcations of period 2 n -tupling on 
the Mandelbrot set. This fact as well explains the 
similarity of the critical Julia set and the Mandel- 
brot set for the original map. 

The sets of unstable cycles of periods 6 — 11 are 
shown in Fig.4 (a-f). One can see that with in- 
crease of n, the full picture of critical set becomes 
more and more distinguishable. The analogy be- 
tween elements of the cycles and Yang-Lee zeros, 
and also between limit of the infinite period of cy- 
cles n and a thermodynamic limit in the theory 
of phase transitions is evident. 
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Figure 4: Sets of all unstable cycles with peri- 
ods n = 6 (a), 7 (b), 8 (c), 9 (d), 10 (e), 11 (f) 
obtained by backward iterations of RG transfor- 
mation (12). 

5 Electrostatical analogy — 
the potential of critical set 
and its properties near the 
transition to chaos 

The phase transition critical point corresponds 
to intersection of the real axis with the critical 
line in the complex plane, formed in the thermo- 
dynamical limit by the set of zeros of partition 
function. This critical point should provide the 
jump of a derivative of the thermodynamic free 
energy (in the case of I-type phase transition) or 
of the second derivative (in the case of the II-type 
phase transition). Let us notice that the distribu- 
tion of free energy depending on temperature can 
be interpreted as the distribution of the electro- 
static potential created in 2D space by the set of 
charges, located at the partition function zeros, 
and in thermodynamic limit - the potential of 
the corresponding limiting distribution of charge 
on the critical Julia set. 

The electrostatic potential for critical Julia 



6 



set can be calculated from the Hubbard and 
Douady |3| formula 



U 



lim ln|<^ n (A)|. 



(16) 



This formula arises from reasons of existence of 
conformal mapping of the Julia set to a circle. 
The factor is not so important. In our case it is 
defined by the third order of the map (12): By 
one iteration, the potential of the point increase 
three times. At infinity, the potential of the Julia 
set coincides with potential of the charged disk 
U = (l/3 n ) ln|4 

We have examined distribution of the electro- 
static potential and its derivatives, both in com- 
plex area, and on the real axis. In Fig. 5 the 
set of several equpotential lines of homogeneously 
charged 2D object limited by Julia set of map 
(12) are shown. This figure characterizes the dis- 
tribution of potential in the complex plane. In 
Fig. 6 the distribution of the electrostatic normal 
component of a field (derivative of the potential) 
along the real axis is shown. The fact of exis- 
tence of the jump of the derivative of the potential 
at the critical point proves to be true. Near the 
other point of crossing of the critical set with the 
real axis the first derivative of potential behaves 
continuously, but has a break providing jump of 
the second derivative. This point has to be in- 
terpreted as the II-type phase transition. As ap- 
pears, the dependence of potential close to the 
Feigenbaum critical point behave asymptotically 
as U(Re\) = (ReA — ReA*) 7 , where the power 
factor 7 ~ 0.63 (see Fig. 7). 

According to these results, we can assert that 
the potential of the Julia set of the RG transfor- 
mation is a certain criterion of ordering of chaotic 
dynamics near to the point of transition to chaos, 
like the free energy is connected with the order pa- 
rameter (magnetization) of thermodynamic sys- 
tem. One of known analogue of the order param- 
eter is Lyapunov exponents [2]. However some- 
times it is more effective to use the electrostatic 
potential for the description of dynamics, for ex- 
ample when the trajectories escape to infinity. 




Figure 5: Julia set of the map (12) and its equpo- 
tential lines, demonstrating distribution of elec- 
trostatic field at the areas of escape to infinity on 
the complex plane. 




Figure 6: Distribution of normal component of 
intensity of the electrostatic field of the critical 
set along the real axis A. 




Figure 7: Distribution of inclination of a double 
logarithmic curve In [/(ReA) vs ln(ReA — ReA*) 
along real axis A near the Feigenbaum critical 
point A*. 
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6 RG analysis of period- 
multiplication bifurca- 
tions cascades 

Let us generalize the approximate renormaliza- 
tion procedure for other bifurcation cascades, dis- 
tinct of that of Feigenbaum. Namely, we consider 
the period-multiplications such as period-tripling, 
period-quadrupling etc. intrinsic to complex an- 
alytic maps. 

To derive the RG transformation of the period- 
tripling (quadrupling) critical point, one must 
apply the original map (J1J) three (four) times 
and represent the result as a Taylor series up to 
quadratic terms in x n : 

x n+3 = /(/(/(*„))) = A- (A -A 2 ) 2 - 

4A(A-A 2 )^ + 0(4), [U) 

x n+ , = f(f(f(f(x n )))) = A - (A - (A - A 2 ) 2 ) 2 + 
8A(A - A 2 )(A - (A - A 2 ) 2 )x 2 + 0{x 4 n ). 

(18) 

From (fTTj) and (|T5j) one can obtain the pa- 
rameter renormalization transformations 
for period tripling (Eq. fl9|) ) and period- 
quadrupling (Eq. PU|) ) respectively: 

A' — > 4A(A - A 2 )(A - (A - A 2 ) 2 ), (19) 

A' -> -8A(A-A 2 )(A-(A-A 2 ) 2 )(A-(A-(A-A 2 ) 2 ) 2 ) 

(20) 

and respective renormalization factors are ex- 
pressed as: 

a = 4A(A-A 2 ), (21) 
« = -8A(A-A 2 )(A-(A-A 2 ) 2 ). (22) 

Julia sets of the transformations (fT9*J) and (|2*nj) 

are shown in figure 8 (a, b). Critical values of 
parameter A, scaling factors, parameter scaling 
constants and critical multipliers are numerically 
calculated. It is evident that the critical indices 
are represented by complex numbers. It is pos- 
sible to compare them with the data of the nu- 
merical RG analysis f2Uj (see Tables 2 and 3). 
For example, let us consider one of the fixed 
points of the transformation ()19j) . which cor- 
responds to the Golberg-Sinai-Khanin point of 
the period-tripling accumulation j^U 122]- Fig- 
ure 9 (a) demonstrates scaling properties of the 
critical set of the transformation (fT9*|) near the 
approximate period-tripling accumulation point. 




Figure 8: Charts of initial values of parameter A of 
RG transformations for different order of bifurca- 
tion cascade N = 3 (a), N = A (b), N = 2*3 (c), 
N = 2 * 2 * 3 (d) 




Figure 9: Scaling of the Julia set of approxi- 
mate RG transformation (a) and Mandelbrot set 
of the map (0J) (b) near the critical point of period- 
tripling bifurcations cascade. 
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Table 2: Rigorous and approximate values of crit- 
ical indexes for the period-tripling accumulation 
point 





Rigorous 


Approximate 


A* 


0.024+0.784i 


0.025+0.792i 


a 


-2.097+2.358i 


-2.317+2.147i 


5 


4.600+8.98H 


7.078-7.624i 




-0.476-1.055i 


-0.493-1. 062i 



Table 3: Rigorous and approximate values of crit- 
ical indexes for the period-quadrupling accumula- 
tion point 





Rigorous 


Approximate 


A* 


-0.310+0.495i 


-0.314+0.493i 


a 


-1.131+3.260i 


-1.238+2.913i 


5 


-0.853-18. HOi 


3.041-13.984i 




0.063-1. 053i 


0.070-1. 060i 



On panel (b) an illustration of the scaling prop- 
erties of the Mandelbrot set near the Golberg- 
Sinai-Khanin point is shown. A small fragment of 
the critical set (fig. 9a) (Mandelbrot set (fig.96)) 
by the multiplication to the approximate (rigor- 
ous) scaling constant 6 transforms to a picture in 
the right column. As the constant 5 is a com- 
plex number, this transformation includes a scale 
change and a rotation. 

Let us generalize the approximate RG analysis 
to the period iV-tupling bifurcations cascades for 
an arbitrary N. By induction, the expression for 
the scaling factor and the parameter renormaliza- 
tion transformation look as follows: 

JV-l 

a = (-2)"" 1 J] /i(0), (23) 

i=l 
N 

A'-H^n^O). (24) 

Obviously, the critical Julia sets of the RG trans- 
formations (|2*4^1 with increasing N approximate 
the fractal properties of the Mandelbrot set more 
and more precisely. For example, one of the 
fixed points of the RG transformation of period- 
quadrupling takes place on the real axes at A* = 
1.396. This point corresponds to the Feigenbaum 
critical point A* = 1.401 and is approximated in 
this renormalization scheme more precisely than 
previous estimate A* = 1.37. 




Figure 10: Julia set of the approximate RG trans- 
formation of bifurcation cascade of infinite order 
N -> oo. 

It is worth noting that for the better approxi- 
mation of the Mandelbrot set by Julia sets of the 
RG transformation, for a large period multiplica- 
tion the factor iV (increasing to infinity) must be 
a composite number, that is iV = 2 • 3 • ... Then 
the considered RG transformation will describe 
a number of different bifurcations cascades (see 
figure 8). 

Let us construct a map, which allows an RG 
transformation of order N + 1 using RG transfor- 
mations of low orders. It yields 

Aam-i = (—2)\ N * 

*(^fc"(^fc) 2 -(^) 2 ), (25) 
Ai = /\(0) = A, A = 1, A_i = oo. 

Dynamics of such a map with N — > oo in- 
cludes description of all period iV-tuplings cas- 
cades. Figure 10 shows the plane of initial values 
of Ai = A. One can see that it precisely corre- 
sponds to the Mandelbrot set. 

7 Conclusion 

In the present paper the complex variable version 
of approximate RG method for period-doubling 
bifurcation cascade is considered. It is shown, 
that the Julia set of renormalization transforma- 
tion of parameter of complex logistic map is the 
approximate version of Mandelbrot set of this 
map. This similarity is explained by the fact, that 
Julia set is a set of every possible unstable cycles 
of RG transformation, and elements of these cy- 
cles correspond to accumulation points of various 
bifurcations of period 2 n multiplication (n — > oo), 
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located on the boundary of Mandelbrot set dense 
everywhere. 

Approximate RG analysis is generalized to the 
case of different cascades of period-multiplication 
bifurcations (for example period-tripling, period- 
quadrupling etc.), which is peculiar for the com- 
plex analytic maps. It is shown that with increase 
of order of the bifurcation cascade the similarity 
between Mandelbrot set and Julia set of renor- 
malization transformation becomes more clear. 

The obtained outcomes are interpreted in a 
view of analogy with theory of phase transitions, 
namely the Yang-Lee theory, based on investiga- 
tion of properties of thermodynamic values de- 
pending on complex temperature. It is necessary 
to note, that the elements of considered unsta- 
ble cycles of RG transformation which is equiv- 
alent to points of Julia set of this transforma- 
tion, in thermodynamic analogy corresponds to 
so-called Yang-Lee zeros, defining in a thermody- 
namic limit the borderline of phase transition. It 
is shown, that at the points of transition to chaos 
the jump of an electrostatic field of a critical Julia 
set is observed. Within the framework of consid- 
ered analogy, the jump of a derivative of a free en- 
ergy at phase transition points takes place. Thus, 
it is possible to define a new criterion of transi- 
tion to chaos. The electrostatic potential can be 
regarded as order parameter for the transition. 

We conclude that complex generalization of ap- 
proximate RG method appears to be useful to ad- 
vance understanding of the critical phenomena at 
threshold of chaos and for development of anal- 
ogy with the theory of phase transition, which 
can give new approaches to investigation of these 
phenomena. 
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